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Section  I 


INTRODUCTION 

It  is  well  recognized  that  the  procedures  employed  during  the  manufacture  of  composite 
materials  have  a  significant  effect  on  the  quedity  and  cost  of  the  finished  product.  Therefore, 
the  procedures  used  during  manufacture  must  carefully  be  selected  for  each  application. 
In  particular,  attention  must  be  paid  to  the  processing  variables,  such  as  the  temperature 
and  pressure,  applied  during  fabrication.  The  processing  variables  can  best  be  chosen  by 
the  use  of  analytic  models.  For  this  reason,  process  models  applicable  to  thermoset  matrix 
composites  have  been  developed  in  recent  years  [1-4].  However,  a  complete  model  of  the 
manufacture  of  thermoplastic  matrix  composites  has  not  yet  been  reported.  Hence,  the 
objective  of  this  investigation  was  to  develop  a  model  simulating  the  major  steps  involved 
in  the  processing  of  semicrystalline  thermoplastic  matrix  composites. 

Processing  of  thermoplastic  matrix  composites  consists  of  three  major  steps: 

1)  Matrix  is  introduced  into  the  space  between  the  fibers  in  each  tow,  “impregnating”  the 
tow. 

2)  Individual  plies  are  “consolidated”  such  that  there  is  good  contact  as  well  as  good  bond 
between  adjacent  phes. 

3)  The  composite  is  heated  and  then  cooled  at  a  rate  which  provides  the  desired  “crys¬ 
tallinity”  in  the  matrix. 

In  this  paper  a  model  is  presented  describing  each  of  the  above  three  steps.  The 
model  is  divided  into  the  three  submodels  designated  as  “impregnation”,  “consolidation”, 
and  “crystallinity”.  Each  of  these  submodels  is  developed  separately.  The  three  submodels 
are  then  combined,  so  as  to  be  applicable  in  situations  where  some  or  all  of  the  phenomena 
described  by  the  different  submodels  occur  simultaneously. 

The  submodels  were  verified  by  data  obtained  with  PEEK  150P  polymer,  with  T300 


1 


fiber,  and  with  APC-2  graphite/PEEK  composite.  It  is  recognized  that  the  grade  of  PEEK 
used  in  APC-2  may  differ  from  PEEK  150P  in  molecular  weight  distribution  and  in  additive 
content.  PEEK  150P  was  chosen  on  the  recommendation  of  Fiberite  Corporation,  an  ICI 
subsidiary.  According  to  the  information  provided  by  Fiberite  Corporation,  150P  is  the 
closest  commercial  grade  of  PEEK  to  the  APC-2  material. 
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Section  II 


IMPREGNATION  SUBMODEL 

The  first  step  in  the  manufacture  of  thermoplastic  composites  is  impregnation.  During 
this  step  a  polymeric  matrix  is  introduced  into  the  fiber  tow.  Impregnation  is  done  either 
by  the  material  supplier  (manufacturer)  during  preparation  of  prepreg  fibers,  plies,  or  lam¬ 
inates,  or  by  the  user  during  “in  situ”  fabrication  of  laminates  and  parts.  Impregnation  is 
accomplished  by  surrounding  the  tow  with  the  matrix  and  heating  the  tow-matrix  system 
to  a  temperature  Tm  at  which  the  matrix  becomes  “soft”,  capable  of  penetrating  into  the 
space  between  the  fibers  (Figure  1).  Our  objective  is  to  determine  a)  the  depth  to  which  the 
matrix  penetrates  in  a  given  time,  and  b)  the  time  required  to  completely  fill  (impregnate) 
the  tow  with  the  matrix. 

Here  we  postulate  that  the  impregnation  process  takes  place  by  the  following  mecha¬ 
nism.  The  fibers  in  the  tow  are  not  perfectly  straight.  Therefore,  the  distance  between  fibers 
varies  along  the  fiber  length  (Figure  2).  At  a  point  where  two  fibers  are  sufficiently  close 
the  matrix  around  one  of  the  fibers  comes  into  contact  with  the  adjacent  fiber,  bridging  the 
gap  between  the  two  fibers.  Driven  by  surface  tension,  the  matrix  then  moves  lengthwise 
along  the  fibers  until  another  bridging  point  is  reached.  The  lengthwise  flow  along  the  fibers 
is  accompanied  by  radial  flow  across  the  fibers.  The  lengthwise  flow  is  modeled  as  laminar 
flow  in  an  annular  channel  of  length  L,  height  /i,  and  inside  radius  r  (Figure  3).  The  radial 
flow  is  modeled  as  flow  through  a  porous  medium. 

By  assuming  a  linear  pressure  drop  along  the  channel,  the  mass  flow  rate  iha  lengthwise 
along  the  channel  may  be  expressed  ais  [5] 

(1) 

where  p  and  p.  are  the  density  and  viscosity  of  the  matrix  at  temperature  Tm,  Xm  is  the 
length  of  the  channel  filled  by  the  matrix,  and  Pco  and  Pc  are  the  pressures  at  the  two  ends 
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DIAMETER:  d 
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MATRIX 


Figure  1.  Illustration  of  the  impregnation  of  a  fiber  tow. 
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of  the  channel  (Figure  3).  The  pressure  Pc  is  related  to  the  pressure  in  the  empty  channel 
Pa  through  the  surface  tension  a.  We  assume  that  the  channel  is  open  to  the  ambient  so 
that  the  pressure  Pa  is  constant.  The  pressure  difference  across  the  meniscus  can  then  be 
written  as  [5] 

2  2 

Pa  — Pc=  (2) 

h  h 

where  9  is  the  wetting  angle  and  a*  =  a  cos  0. 

The  mass  flow  rate  in  the  radial  direction  through  an  element  of  length  dx  is 


dthr  =  pur{2irr)dx 


(3) 


Ur  is  the  average  radial  velocity  given  by  Darcy’s  law  [5] 

p  dr 


(4) 


Kp  is  the  permeability  perpendicular  to  the  tow,  and  can  be  approximated  by  the  expression 
given  in  Table  1.  Equations  (3)  and  (4)  yield 

2tpK 


TWr  = 


/iln 


^  {P-Pa)dx 
To  Jo 


(5) 


where  P  is  the  pressure  at  a  distance  x  along  the  channel.  Pa  is  the  matrix  pressure  at  the 
outside  radius  ro  of  the  tow,  and  is  assumed  to  be  the  same  as  the  pressure  which  exists 
inside  the  empty  portion  of  the  channel. 

Conservation  of  mass  requires  that  the  lengthwise  and  axial  mass  flow  rates  be  equal 


ma  —  rrir 


(6) 


By  combining  Equations  (1)  through  (6),  after  lengthy  but  straightforward  algebraic  manip¬ 
ulations,  we  obtain 

dxm  _  2<t*  Kp _ _ 

dt 

Initially  {t  =  0)  there  is  a  small  amount  of  matrix  on  the  bottom  of  the  channel.  We  take 
the  length  of  this  amount  to  be  equal  to  the  channel  height  h.  The  time  required  to  fill 
the  channel  with  matrix  {xm  =  L)  is  denoted  by  ti.  Thus  we  integrate  Equation  (7)  from 
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Table  1 

INPUT  PARAMETERS  REQUIRED  FOR  THE  IMPREGNATION  SUBMODEL 


1)  Initial  tow  radius,  Tq 

2)  Fiber  diameter,  d 

3)  Average  distance  between  two  fibers  (channel  height),  h 

4)  Channel  length,  L 

5)  Permeability  of  the  tow,(“)  Kp 

6)  Matrix  viscosity,  fi 

7)  Surface  tension, a* 

8)  Applied  temperature,  Tm 


a)  The  permeability  may  be  approximated  by  [6] 


where 


b)  In  the  temperature  range  of  interest  the  surface  may  be  expressed  as 


<T*  —  <T*e°‘ 


T-T, 


where  cr*  and  a  axe  constants  which  depend  on  the  matrix,  fiber,  and  interface, 
but  are  independent  of  temperature.  To  is  an  arbitrary  reference  temperature. 
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xm  =  h  to  Xm  =  L  and  from  t  =  0  to  t  =  tL.  This  integration  yields  the  following  expression 
for  the  time  required  to  fill  a  channel  located  at  radius  r 


tl  =  G  1  —  /rln 


G  and  /  are  constants  defined  as 


(31 


Zuh  [{L/hf  - 1] 


2a* 

^  h  \n{L/h) 

'-3/4[(W-l] 

It  is  interesting  to  note  that  Equation  (8)  is  also  obtained  if  the  pressure  along  the  fibers  is 
taken  to  be  constant  instead  of  varying  linearly  with  position,  as  is  assumed  in  the  above 
analysis.  However,  for  constant  pressure  the  value  of  G  would  be  increased  while  the  value 
of  /  would  be  decreased  by  a  factor  of  two. 

The  time  required  to  fill  the  tow  to  a  depth  of  r  (i.e.,  the  time  required  to  fill  every 
channel  between  Vg  and  r  (Figure  1))  is 

( = (11) 
i=l 

where  Nr  is  the  number  of  “layers”  between  r  =  r<,  and  r  =  r.  The  thickness  of  each  layer  is 

Ar  =  d  +  h  (12) 

where  d  is  the  fiber  diameter.  Thus,  the  number  of  layers  is  Nr  —  {f'o  r)/Ar.  We  can 
replace  the  summation  in  Equation  (11)  by  the  integral 


Substitution  of  Equation  (8)  into  Equation  (13)  and  integration  yields 


t*  =  {l-  r*)  +  Jr*^]nr*  +  -(1  -  r*^) 


t*  and  r*  are  dimensionless  time  and  position,  respectively,  and  J  is  a  dimensionless  constant 

_  ^^(0  -  2<T*(d/h  + 1)  ^ 

Gro  3/i  [{L/h)^  -  l] 


T  _  1^0  _ 


r*  =  r/ro 


hro  ]ii{L/h) 
f,  [{L/kf  - 1' 
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Equation  (14)  provides  the  time  required  to  impregnate  the  tow  with  matrix  to  a  depth  of 
(ro  —  r).  It  is  convenient  to  express  the  extent  of  impregnation  by  a  parameter  referred  to 
as  the  degree  of  impregnation 


^  _  area  impregnated  _  7r(r^  ~  ^ 

■^imp  —  :  :  T~!  —  - ^ -  =  1  ~  1  — 


total  tow  area 


7rrf 


In  terms  of  the  degree  of  impregnation  Equation  (14)  is 


(18) 


t*  =  (l-  v/1  -  Amp)  +  ^(1  -  i^imp)  lny/l-Di^j,  +  -D, 


2  -  imp 


The  time  required  to  fill  the  entire  tow  is  (r  =  0,  and  Amp  =  1) 


Cp=l+2 

or,  in  terms  of  dimensional  quantities,  the  time  to  complete  the  impregnation  is 

hro]n(L/h) 


*'imp 


_  3firo  [(L/h)^  ~  ^1  fi 
2cr*(d/h 1)  \ 


+ 


12Kj,[(L/h)2 


-Hi 


(19) 


(20) 


(21) 


Impregnation — Method  of  Solution 


The  degree  of  impregnation  as  a  function  of  time  may  be  calculated  by  Equation  (19). 
The  time  required  to  complete  the  impregnation  may  be  calculated  by  Equation  (20)  or  (21). 
The  calculations  can  be  facilitated  by  representing  Equations  (19)  and  (20)  in  graphical  form. 
Graphical  representation  of  these  equations  are  given  in  Figures  4  and  5.  Since  these  figures 
are  in  terms  of  dimensionless  parameters  they  are  not  restricted  to  specific  materials  or 
specific  tows  but  are  applicable  to  different  types  of  materials  and  different  types  of  tows. 
The  actual  degree  of  impregnation  and  the  actual  impregnation  time  can  be  calculated 
from  these  figures  once  the  parameters  applicable  to  the  given  problem  are  known  (see 
Table  1).  Of  the  parameters  required  to  calculate  Dimp  and  fimp,  the  tow  radius  To,  the  fiber 
diameter  d,  the  distance  between  the  fibers  (channel  height)  h,  and  the  channel  length  L 
can  be  determined  from  visual  observation  of  the  tow.  The  resin  viscosity  /i  as  a  function 
of  temperature  can  be  measured  by  a  viscometer.  The  permeability  Kp  can  be  estimated  by 
the  expression  given  in  Table  1.  The  surface  tension  <7*  must  be  determined  by  matching 
the  model  to  data  obtained  by  impregnation  tests.  A  procedure  for  these  tests  is  described 
in  the  next  subsection. 
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Figure  4.  Degree  of  impregnation  as  a  function  of  time  (dimensionless,  see  Equation  (19)). 


ONLESS  PARAMETER,  J 


regnate  a  tow  characterized  by  the  dimensionless 

(20)). 


The  paxeimeters  required  to  calculate  the  degree  of  impregnation  and  the  impregnation 
time  were  measured  for  T300/PEEK  150P  (see  next  subsection).  The  measured  values  of  the 
parameters  are  hsted  in  Appendix  A.  With  these  parameters,  the  degree  of  impregnation  as 
a  function  of  time  and  the  time  of  impregnation  as  function  of  tow  size  were  calculated  and 
plotted  for  impregnation  occurring  at  698®F  (370®C)  (Figures  6  and  7).  The  time  required 
to  complete  the  impregnation  at  different  values  of  were  also  plotted  (Figure  8). 

Finally,  it  is  noted  that  Equations  (19)  through  (21)  were  incorporated  into  a  computer 
code  which  combines  the  calculations  of  all  three  submodels  (impregnation,  consolidation, 
crystalhnity).  This  code  is  discussed  in  Section  V. 

Impregnation — Experimental  Verification 

Tests  were  performed  to  evaluate  the  validity  of  the  impregnation  submodel,  to  estabhsh 
the  procedure  for  determining  some  of  the  parameters  required  by  this  submodel,  and  to 
measure  the  parameters  for  one  fiber-matrix  system.  In  these  tests  Union  Carbide  T300 
graphite  fiber  tows  (6K)  were  impregnated  w'ith  PEEK  150P  polymer.  Prior  to  the  tests  the 
sizing  was  removed  from  the  fibers  by  the  following  procedure.  The  fibers  were  immersed  in 
methylethylketone  (MEK)  for  24  hours,  then  were  washed  with  fresh  MEK  in  an  ultrasonic 
cleaner  for  30  minutes.  This  process  was  repeated.  After  the  second  treatment  with  fresh 
MEK  the  fibers  were  dried  in  a  vacuum  oven  at  180®F  (82.2°C)  for  two  hours. 

A  fiber  tow,  treated  in  the  above  manner,  was  placed  in  an  aluminum  mold  (Figure  9). 
The  two  ends  of  the  tow  were  clamped  so  that  the  tow  was  suspended  above  the  bottom  of 
the  mold.  The  mold  was  filled  with  PEEK  150P  powder  and  placed  in  an  oven  kept  at  698'’F 
(370°C).  The  temperature  of  the  matrix  was  monitored  by  a  thermocouple.  After  a  preset 
period  of  time  the  mold  was  removed  from  the  oven  and  cooled  to  room  temperature.  The 
impregnated  tow  was  removed  from  the  mold,  was  cut  perpendicular  to  the  fibers,  and  the  tow 
cross  section  was  examined  with  an  optical  microscope.  The  number  of  impregnated  fibers 
in  the  tow  was  counted,  and  the  degree  of  impregnation  was  calculated  by  the  expression 


•^imp  — 


number  of  impregnated  fibers 
total  number  of  fibers 


(22) 
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TIME,  t  (hours) 


Figure  6. 


Degree  of  impregnation  as  a  function  of  time  as  calculated  by  the  model  ( 
Equation  (19)). 
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TIME  TO  COMPLETE  IMPREGNATION,  t  imp 


seconds  hours 


Figure  8.  Time  required  to  impregnate  a  tow  cis  calculated  by  the  model  (see  Equa¬ 
tion  (21)). 
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Figure  9.  Schematic  of  the  mold  iised  in  the  impregnation  tests. 


This  process  was  repeated  with  tows  having  been  in  the  oven  for  different  lengths  of  time. 
In  this  manner  the  degree  of  impregnation  as  a  function  of  time  was  determined.  The  data 
are  shown  in  Figure  10. 

The  degree  of  impregnation  was  also  calculated  by  the  impregnation  submodel  (Equa¬ 
tion  19)  for  the  conditions  of  the  tests.  The  geometric  parameters  required  in  the  calculations 
(channel  height  h,  tow  radius  r^,  and  fiber  diameter  d)  were  measured  directly  from  pho¬ 
tomicrographs  of  the  tow  cross  section  perpendicular  to  the  fibers. 

The  average  channel  length  L  was  determined  from  photomicrographs  of  tow  sections 
parallel  to  the  fibers.  This  was  done  in  the  following  manner.  First,  the  angles  between 
different  pairs  of  fibers  were  measured,  and  from  these  measurements  an  average  angle  Q 
was  computed.  Second,  the  average  channel  length  was  calculated  by 

^  _  channel  height  h 
tan  ft  tan  ft 

The  viscosity  of  P EEK  150P  was  measured  with  a  Rheometrics  parallel  plate  type  viscometer. 
The  value  of  the  surface  tension  a*  was  determined  by  fitting  the  model  to  data.  The  values 
of  the  material  properties  are  given  in  Appendix  A. 

The  degree  of  impregnation  calculated  by  the  model  is  included  in  Figure  10.  There 
is  reasonable  agreement  between  the  measured  and  calculated  degree  of  impregnation,  and 
this  lends  support  to  the  validity  of  the  model. 
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Figure  10.  Degree  of  impregnation  as  a  function  of  time.  Comparisons  of  data  and  the  results  of  the  model. 


Section  III 


CONSOLIDATION  SUBMODEL 

When  thermoplastic  composites  are  being  processed  individual  plies  consolidate  into  a 
laminate  by  bonding  at  the  interfaces  [7].  This  bonding  consists  of  two  phenomena.  First, 
two  adjacent  ply  surfaces  coalesce  and  come  into  “intimate  contact”.  Second,  bond  forms 
at  the  ply  interface  by  a  process  called  autohesion.  In  the  following,  models  are  presented 
which  describe  the  effects  of  the  processing  variables  (temperature,  pressure,  time)  on  the 
intimate  contact  and  autohesion  processes.  These  processes  occur  simultaneously.  However, 
for  the  sake  of  clarity,  the  models  of  these  processes  are  developed  separately. 

Intimate  Contact — Model 

Laminates  are  formed  by  “laying  up”  plies.  Since  the  ply  surfaces  axe  uneven,  spatial 
gaps  exist  between  the  plies  prior  to  the  application  of  heat  and  pressure.  When  heated,  most 
thermoset  matrix  composites,  including  epoxy,  have  sufficiently  low  viscosities  and  wetting 
abilities  to  coalesce  the  ply  surfaces.  For  thermoplastics,  even  at  elevated  temperatures,  the 
viscosity  is  too  high  to  produce  the  desired  degree  of  flow  along  the  interface.  Therefore, 
thermoplastic  matrix  composites  must  actually  be  deformed  to  produce  intimate  contact 
between  adjacent  surfaces. 

Following  Dara  and  Loos  [7]  we  represent  the  irregular  ply  surface  by  a  surface  con¬ 
sisting  of  a  series  of  rectangles  (Figure  11).  However,  while  in  their  model  Dara  and  Loos 
utilized  rectangles  of  different  sizes,  we  take  the  rectangles  to  be  identical.  Due  to  the  applied 
force  (or  pressure)  the  rectangular  elements  spread  along  the  interface.  We  are  interested  in 
the  extent  of  this  spread  (degree  of  intimate  contact)  as  a  function  of  applied  temperature, 
force  (pressure),  and  time.  We  also  wish  to  calculate  the  time  required  to  reach  “complete” 
intimate  contact. 
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F=  (Pressure)  x( Area) 


I 


F 


Figure  11.  Illustration  of  the  idealized  interface  used  in  the  intimate  contact  model 
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With  reference  to  Figure  12,  the  degree  of  intimate  contact  is  defined  as 


b 

Wo  +  bo 


(24) 


where  bo  and  6  are  the  initial  {t  <  0)  and  instantaneous  (at  time  t)  widths  of  each  rectangular 
element,  respectively,  and  Wo  is  the  initial  distance  between  two  adjacent  elements.  During 
processing  the  volume  of  each  element  remains  constant 


Vo  —  Oiobo  — 


(25) 


where  ao  and  a  are  the  initial  and  instantaneous  heights  of  each  rectangular  element.  Equa¬ 
tions  (24)  and  (25)  yield  the  following  expression  for  the  degree  of  intimate  contact 


D\c  = 


0,0  !  o, 


1  -b  Wojbo 


(26) 


To  proceed  with  the  model  we  apply  the  law  of  conservation  of  mass  to  a  control  volume  of 
width  dy.  Referring  to  Figure  12  we  write 


duy  do  _ 

o — -  H - =  0 

dy  dt 


(27) 


y  is  the  coordinate  along  the  interface  and  t  is  time.  By  assuming  that  the  flow  is  laminar, 
we  write  the  average  velocity  Uy  as  [5] 


-  dP 

“  12/im/  dy 


(28) 


Here  Hmf  is  the  viscosity  of  the  fiber- matrix  mixture.  In  the  space  between  two  adjacent 
elements  the  pressure  is  Pe.  The  edge  of  the  element  [y  =  b/2),  moves  with  a  speed  of  dbfdt. 
Thus,  the  boundary  conditions  corresponding  to  Equation  (28)  are 


P  =  Pe 


and 


db 


at  y  =  6/2  for  t  >  0 


(29) 


By  combining  Equations  (27),  (28),  and  (29),  after  algebraic  manipulations,  we  obtain 


P  —  P  —  ^ 

a3  dt 


(30) 
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Figure  12.  Top:  Rectangular  elements  representing  the  uneven  surface  at  time  t  —  0. 

Bottom:  Illustration  of  one  element  at  time  t,  and  the  control  volume  used  in 
calculating  mass  flows. 


The  force  applied  to  the  entire  ply  of  length  A  and  width  B  is  F  (Figure  11).  Correspondingly, 
the  force  applied  per  unit  length  to  one  element  is 


(31) 


where  n  -  Bf  (bo+Wo)  is  the  number  of  elements  in  a  ply  of  width  B.  We  write  Equation  (31) 


as 


f  -  L  +  ^o)  _  p  (h  ,  A 


(32) 


where  P^pp  is  the  applied  gauge  pressure,  and  P*  is  the  ambient  pressure  which  is  taken  to  be 
the  same  as  the  pressure  in  the  space  between  the  rectangular  elements.  The  force  applied 
to  an  element  must  be  balanced  by  the  pressure  inside  the  element 

rt/2 

I 

-6/2 

Combination  of  Equations  (25),  (30),  (32),  and  (33)  and  integration  yields 


rb/2 

/=/  (P-Pe)dy 
J-b/2 


(33) 


Oo 

a 


i  +  1 


Kj  \bj 


2  -|l/5 
t 


(34) 


Equations  (34)  and  (26)  provide  the  following  expression  for  the  degree  of  intimate 
contact 

1/5 

(35) 


Die  = 


1 


1  +  ^ 


j  ^  ^app 


(■*5)  te) 


2  1 
t 


Complete  intimate  contact  is  achieved  when  Die  becomes  unity.  The  height  a 
spending  to  this  condition  is  (see  Equation  26) 


corre- 


Substitution  of  Equation  (36)  into  Equation  (35)  results  in  the  following  expression  for  the 
time  required  to  achieve  complete  intimate  contact 

2 


t'lr  — 


5P 


app 


l+Wo/bo  \ao  J  \  bo  J 


(37) 
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Intimate  Contact — Method  of  Solution 


The  degree  of  intimate  contact  and  the  time  required  to  reach  complete  intimate  contact 
can  be  calculated  by  Equations  (35)  and  (37),  respectively.  The  parameters  required  for  the 
calculations  are  summarized  in  Table  2.  The  geometric  parameters  WolK  and  aojho  can  be 
measured  from  photomicrographs  of  the  cross  section  of  an  uncompacted  ply.  The  viscosity 
Hmf  can  be  obtained  by  matching  the  model  to  degree  of  intimate  contact  versus  time  data. 
A  method  for  generating  such  data  is  given  in  the  next  subsection. 

Values  of  WolK,  OolK,  and  fXmf  for  APC-2  are  included  in  Appendix  A.  With  these 
values  the  degree  of  intimate  contact  Ac  and  the  time  to  achieve  complete  intimate  contact 
tic  were  calculated  (solid  lines  in  Figures  13  and  14).  Similar  plots  can  be  generated  for 
other  types  of  materials  once  the  values  of  the  relevant  parameters  are  known.  Methods  for 
determining  these  parameters  are  given  in  the  next  subsection. 

Equations  (35)  and  (37)  were  also  incorporated  into  a  computer  code  which  simulates 
the  entire  manufacturing  process  (Section  V). 

Intimate  Contact — Experimental  Verification 

The  intimate  contact  model  was  evaluated  by  the  following  test  procedure.  The 
aluminum- steel  mold  shown  in  Figure  15  was  placed  into  a  hot  press  and  was  preheated 
to  a  given  temperature.  One  ply  of  APC-2  was  then  placed  between  the  top  and  bottom 
aluminum  plates  of  the  mold  and  a  pressure  was  applied.  After  a  certain  amount  of  time  the 
pressure  was  released  and  the  mold  was  cooled  to  room  temperature.  The  ply  was  then  taken 
out  of  the  mold  and  its  surfaces  were  inspected.  The  areas  of  the  uncompressed  surfaces 
were  measured,  and  from  these  measurements  the  degree  of  intimate  contact  was  determined 

^  (total  surface  area)  —  (uncompressed  area) 
total  surface  area 

These  tests  were  performed  at  40  psig  and  662,  680,  and  698‘’F  (350,  360,  and  370°C),  and 
at  40,  96,  and  227  psig  at  662°F  (350°C).  The  results  of  these  tests  are  included  in  Figure  13. 
There  is  good  agreement  between  the  results  of  the  model  and  the  data  over  the  range  of 
pressures  and  temperatures  used  in  the  tests. 
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Table  2 

INPUT  PARAMETERS  REQUIRED  FOR  THE  CONSOLIDATION  SUBMODEL 


Intimate  Contact 

1)  Length  of  the  plate,  A 

2)  Width  of  the  plate,  B 

3)  Geometric  ratios  aoiwo  and  bg/wo  (see  Figure  12) 

4)  Matrix  fiber  viscosity,  Hmf 

5)  Force  (or  pressure)  applied  to  the  plate,  F  or  ^app 

6)  Applied  temperature,  Tm 

Autohesion 

1)  Proportionality  constant,  k  (see  Equation  41) 

2)  Applied  temperature,  Tm 
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DEGREE  OF  INTIMATE  CONTACT,  D: 


350  °C  (662  °F) 


40psi  (276KPa) 


0  200  400  200  400 


TIME,  t  (sec) 


Figure  13.  Degree  of  intimate  contact  versus  contact  time  as  a  function  of  applied  tem¬ 
perature  and  pressure.  Comparison  of  the  results  of  the  model  with  data. 
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UNITS  IN  INCHES 


Figure  15.  Mold  used  in  the  intimate  contact  tests. 


Autohesion — Model 


Once  two  adjacent  interfaces  come  into  intimate  contact,  a  bonding  process  between 
the  interfaces  starts.  In  the  case  of  two  similar  thermoplastic  interfaces  the  bond  is  formed 
mostly  by  autohesion.  During  autohesion  segments  of  the  chain  like  molecules  diffuse  across 
the  interface  (Figure  16).  The  extent  of  the  molecular  diffusion  and,  hence,  the  bond  strength 
increase  with  time. 


A  convenient  way  of  characterizing  the  extent  of  autohesion  is  through  the  bond 
strength.  Thus  a  degree  of  autohesion  may  be  defined  as  [7] 


= 


(39) 


where  S  is  the  bond  strength  at  time  t  and  Soo  is  the  ultimate  bond  strength,  i.e.,  the 
strength  of  a  completely  bonded  interface. 


We  approximate  the  autohesion  by  the  following  expression  proposed  by  previous  in¬ 
vestigators  (e.g.,  see  References  7  through  10)  for  amorphous  polymers 


D 


au 


(40) 


where  ta  is  the  time  elapsed  from  the  start  of  the  autohesion  process  (i.e.,  the  time  when  the 
interfaces  come  into  intimate  contact),  k  is  a  constant  which  is  related  to  the  temperature 
through  the  Arrhenius  relation 

K  =  Koexp(-E/RT)  (41) 

Kq  is  a  constant,  E  is  the  activation  energy  and  R  is  the  universal  gas  constant.  It  is 
noted  that  the  crystallinity  of  the  material  resulting  from  cooling  may  affect  the  accuracy  of 
Equation  40. 

The  time  required  to  complete  the  autohesion  process  is  discussed  in  the  next  section. 
First,  experimental  verification  of  the  model  is  presented. 

Autohesion — Experimental  Verification 

The  autohesion  model  was  verified  by  Dara  and  Loos  [7]  for  UDEL  PI 700  polysulfone 
film.  Dara  and  Loos  measured  the  bond  strength  between  two  films  brought  into  contact  for 
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t  =  0 
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Figure  16.  Illustration  of  the  autohesion  process. 


a  period  of  time  at  temperatures  above  the  glass  transition  temperature.  The  semicrystalline 
PEEK  150  polymer  must  be  processed  above  the  melting  temperature  644®F  (340°C).  At 
these  temperatures  the  viscosity  of  PEEK  150P  is  very  low  and,  hence,  the  method  of  Dara 
and  Loos  cannot  be  applied  to  PEEK  150P.  The  method  could  be  applied  to  APC-2.  The 
tests  were  performed  as  follows. 

Two  cross  ply  ([O/QOjj')  APC-2  composite  laminates  were  placed  between  two  steel 
sample  holders  (Figure  17).  A  0.001  inch  thick  aluminum  foil,  coated  with  Frekote  44  was 
placed  between  the  two  composite  laminates.  There  was  a  0.16  inch  diameter  hole  in  the 
center  of  the  aluminum  foil. 

The  assembly  described  above  was  placed  in  a  hot  press  maintained  at  either  698, 
716,  or  734°F  (370°C,  380°C,  or  390°C).  A  light  load  (18  Ibf)  was  applied  to  assure  contact 
between  the  laminates.  After  a  certain  amount  of  time  the  assembly  was  removed  from  the 
press  and  was  quenched  in  cold  water..  The  tensile  strength  of  the  bond  was  then  measured 
by  placing  the  assembly  into  a  mechanical  tester  and  by  applying  a  tensile  load  to  the  steel 
sample  holders.  (The  APC-2  bonded  to  the  steel  without  glue.)  The  degree  of  autohesion 
was  calculated  by 

r  i  L  •  bond  strength  at  time  t 

Degree  of  autohesion  - - - -  (42) 

max  bond  strength  (at  t  — »•  oo)  ^  ' 

The  data  thus  obtained  are  given  in  Figure  18.  The  model  (Dau  ~  is  represented  by 
the  solid  lines  in  this  figure.  There  is  reasonable  agreement  between  the  data  and  the  model. 

The  tests  were  repeated  at  different  pressures.  As  seen  from  the  data  in  Figure  19, 
within  the  pressure  range  tested,  the  applied  pressure  does  not  have  a  marked  effect  on  the 
degree  of  autohesion. 

Finally  it  is  noted  that  the  constant  k  can  be  obtained  by  fitting  Equation  40  to  the 
data.  An  expression  of  «  obtained  in  this  manner  is  given  in  Appendix  A. 

Degree  of  Bond 

Once  intimate  contact  is  established  at  a  point  along  the  interface  autohesion,  and  thus 
the  bonding  process,  starts.  The  degree  of  bonding  £)j  at  time  t  can  be  calculated  by  the 
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Figure  17.  Illustration  of  the  test  used  to  measure  autohesion. 
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Figure  18.  Degree  of  autohesion  as  a  function  of  time  for  different  temperatures 
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Figure  19.  Autohesion  as  a  function  of  pressure. 


expression 


t  “  [(^ic)at  AI  ol  ^ 

+  l(Ac)at  2At  (^i.)a.A.][»(*-2A()‘'‘]  + 


(43) 


or 


t/At 


(fl»)  =  Yi  -  (Ac),i-„A,]  [«(<  - 

1=1 


(44) 


where  At  is  an  arbitrary  time  step.  is  the  degree  of  intimate  contact  given  by  Equa¬ 
tion  (35). 


Bonding  is  complete  when  Dj,  becomes  unity.  Thus  the  time  required  to  complete  the 
bonding  is  to  be  calculated  by 

ti/At 

1  =  ^  [(Ac)iAi  ~  (Ac)(i_i)A<]  ~  (45) 

i=l 

The  degree  of  bonding  Dfr  as  function  of  time,  temperature,  and  pressure,  and  the  bonding 
time  tj,  as  function  of  temperature  and  pressure  were  calculated  for  APC-2.  The  results 
are  given  in  Figures  20  and  21.  From  these  figures  the  degree  of  bonding  and  the  time 
required  to  achieve  full  bonding  can  readily  be  estimated.  To  facilitate  the  calculations 
for  conditions  and  materials  not  included  in  these  figures.  Equations  (44)  and  (45)  were 
also  incorporated  into  the  computer  code  developed  for  simulating  the  entire  manufacturing 
process  (Section  V). 
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Figure  20.  Degree  of  bonding  as  a  function  of  time  for  different  temperatures. 
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PROCESSING  TEMPERATURE 


Section  IV 


CRYSTALLINITY  SUBMODEL 

For  semicrystalline  thermoplastic  matrix  composites  the  degree  of  crystallinity  affects 
significantly  the  mechanical  properties  of  the  composite.  The  desired  crystallinity  is  achieved 
during  the  processing  by  cooling  the  composite  from  the  melt  temperature  Tm  at  an  appro¬ 
priate  coohng  rate.  The  objective  here  is  to  establish  a  model  which  relates  the  cooling  rate 
applied  during  processing  to  the  crystallinity  of  the  material.  For  clarity,  the  model  is  devel¬ 
oped  for  a  flat  plate  in  which  the  temperature  varies  only  across  the  plate  but  not  along  the 
plate  (one  dimensional  problem,  Figure  22).  The  model  could  readily  be  extended  to  more 
complex  geometries.  In  applying  the  model  it  must  be  borne  in  mind  that  the  model  was 
developed  to  aid  in  the  selection  of  the  processing  variables.  No  attempt  was  made  to  study 
the  detailed  molecular  processes  and  the  morphologies  involved  in  processing  thermoplastics. 

During  cooling  the  crystallinity  of  the  polymer  changes.  The  instantaneous  degree 
of  crystallinity  depends  on  the  temperature  and  on  the  rate  of  change  of  the  temperature. 
Therefore,  to  determine  the  degree  of  crystallinity  as  a  function  of  position  and  time,  the 
temperature  distribution  inside  the  laminate  must  be  known  at  all  times.  For  a  flat  plate 
the  temperature  distribution  can  be  calculated  by  the  following  form  of  the  conservation  of 
energy  [11] 


where  t  is  time,  z  is  the  coordinate  normal  to  the  plate,  T  is  the  temperature,  p  is  the  density, 
C  is  the  specific  heat,  and  k  is  the  thermal  conductivity  of  the  composite.  Expressions  for 
estimating  the  latter  three  parameters  are  given  in  Table  3. 

The  second  term  on  the  right-hand  side  of  Equation  (46)  represents  the  heat  generated 
due  to  crystallization.  In  this  term  c  is  the  crystallinity  of  the  matrix,  rrim  is  the  matrix 
mass  fraction,  and  is  the  theoretical  ultimate  heat  of  crystallization  of  the  polymer  at 
100  percent  crystallinity. 
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Figure  22.  Geometry  used  in  the  crystallinity  submodel. 
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Table  3 

INPUT  PARAMETERS  REQUIRED  FOR  THE  CRYSTALLINITY  SUBMODEL 
Geometry 

1)  Number  of  plies  in  the  plate,  N 

2)  Thickness  of  one  ply,  $o 

Polymer^®) 

3)  Mass  fraction,  rum 

4)  Density,  pm 

5)  Specific  heat.  Cm 

6)  Thermal  conductivity,  km 

7)  Theoretical  ultimate  heat  of  crystallization, 

8)  Relationship  between  crystallinity,  cooling  rate,  and  temperature 

Fiber(“) 

9)  Density,  pf 

10)  Specific  heat,  Cf 

11)  Thermal  conductivity,  kf 

Mechanical  Properties 

12)  Mechanical  properties  as  function  of  crystallinity  (optional) 

Processing  Variables 

13)  Applied  temperature  or  heat  flux  {Ta  and  Tg  or  and  qs) 


a)  The  density,  specific  heat,  and  thermal  conductivity  of  the  composite  may  be  approxi¬ 
mated  by  the  expressions  [12,  13] 

P  ~  '^mPm  "h  (I  ^m)  Pf 

C  =  irimCm  +  (1  -  rrim)  Cf 

1 


k  = 


(1  -  ^/T^)  + 


\/l/(l  -Vm)  +  ikm/kf  -  1)J 

where  the  subscripts  m  and  /  refer  to  the  matrix  aind  fiber  respectively.  Vm  is  the 
matrix  volume  fraction 

-1 


Vm  = 


l  +  (P: 
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Equation  (46)  contains  two  dependent  variables,  the  temperature  T  and  the  crys¬ 
tallinity  c.  Therefore,  to  proceed  with  the  solution  an  additional  expression  is  needed  relat¬ 
ing  the  crystallinity  to  the  temperature.  This  expression  is  established  by  utilizing  the  fact 
that  for  a  given  material,  the  rate  of  degree  of  crystallinity  depends  on  the  cooling  rate  and 
on  the  instantaneous  temperature.  This  dependence  can  be  expressed  symbolically  as 


^ 

dt  ^  \dt 


,t) 


(47) 


The  function  5  is  a  material  property  which  must  be  obtained  experimentally.  A  procedure 
for  determining  this  function  is  described  in  Appendix  B. 


Solutions  to  Equations  (46)  and  (47)  require  that  the  initial  and  boundary  conditions  be 
specified.  Initially,  (prior  to  the  start  of  cooling)  the  composite  is  at  the  uniform  temperature 
r,  having  zero  crystallinity.  Thus,  the  initial  conditions  corresponding  to  Equation  (46)  are 


T  =  Ti 

c  =  0 


t  <  0 
0<  z<s 


(48) 


where  s  is  the  thickness  of  the  plate.  During  processing  either  the  surface  temperatures  (T^ 
and  Tb)  or  the  heat  fluxes  {q^  and  qs)  on  the  two  surfaces  of  the  plate  must  be  specified 
(Figure  22).  In  terms  of  the  surface  temperatures  the  boundary  conditions  are 


T  =  Ta  z  =  0  t>0 
T==Tb  z  =  s  t>0 


(49) 


In  terms  of  the  heat  fluxes  the  boundary  conditions  are 

q  =  qA  z  =  0  t  >0 
q  =  qs  z  =  s  t>0 

Equations  (46)  through  (50)  define  the  problem  under  consideration. 


(50) 


Crystallinity — Method  of  Solution 

The  temperature  and  the  crystallinity  as  functions  of  position  and  time  can  be  calcu¬ 
lated  by  Equations  (46)  through  (50).  The  solutions  to  these  equations  require  numerical 
procedures.  Therefore,  a  finite  element  algorithm  was  developed  for  performing  the  calcu¬ 
lations.  This  algorithm  was  incorporated  into  a  computer  code  which  simulates  the  entire 
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manufacturing  process  (Section  V).  The  input  parameters  required  for  the  calculations  are 
summarized  in  Table  3. 

To  illustrate  the  type  of  information  which  can  be  generated  by  the  crystallinity  sub¬ 
model,  temperature  and  crystallinity  distributions  as  functions  of  time,  as  well  as  the  average 
crystallinity  as  a  function  of  cooling  rate  were  calculated  for  APC-2  composites.  The  material 
properties  used  in  the  calculations  are  listed  in  the  Appendices.  The  results  are  presented 
in  Figures  23  through  25.  In  these  figures  the  cooling  rates  indicated  axe  the  cooling  rates 
on  the  plates’  surfaces.  It  is  interesting  that,  according  to  the  calculations,  even  for  the 
comparatively  thick  (s  =  0.2  inch)  plates,  the  crystallinity  is  quite  uniform  across  the  plate 
for  coohng  rates  up  to  1000°F/min.  Only  at  the  extremely  high  cooling  rate  of  10,000°F/min 
is  there  a  small  variation  in  crystallinity  inside  the  plate.  This  is  reassuring  because  uniform 
crystalhnity  implies  uniform  material  properties. 

Finally,  it  is  observed  that,  as  expected,  the  average  crystallinity  of  the  plate  changes 
with  cooling  rate  (Figure  25).  At  cooling  rates  below  lOO^F/min  (i.e.,  at  the  cooling  rates 
often  used  in  practice)  the  average  crystallinity  is  nearly  independent  of  the  plate  thickness. 

By  generating  results  such  as  shown  here,  the  cooling  rates  (and  the  corresponding 
surface  temperatures  or  surface  heat  fluxes)  can  be  selected  which  result  in  a  composite  of 
the  desired  crystallinity. 

Crystallinity — Experimental  Verification 

The  validity  of  the  crystallinity  submodel  described  above  was  assessed  by  two  types 
of  tests.  First,  the  crystallinities  of  PEEK  150P  polymer  specimens  cooled  at  different  rates 
were  measured.  Second,  during  cooling  the  temperatures  inside  APC-2  composites  were 
recorded.  The  data  thus  obtained  were  compared  to  the  results  of  the  model. 

The  apparatus  and  procedures  used  in  the  tests  were  described  in  References  14-16. 
Hence,  details  of  the  experiments  are  not  given  here.  Essentially,  in  the  tests,  the  spec¬ 
imens  were  heated  to  720°F  (382'’C)  and  then  cooled  to  room  temperature.  The  surface 
temperatures  of  the  specimens  Ta  were  measured.  An  average  cooling  rate  was  computed  by 

_  (75)at  720® F  ”  (^5)80®F 

^  (At) 
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AVERAGE  CRYSTALLINITY,  c  (weight  percent) 


COOLING  RATE,  Q  (°F/min) 


Figure  25.  Average  crystallinity  as  a  function  of  surface  cooling  rate  for  an  APC-2  com¬ 
posite.  Results  of  the  model. 
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where  Af  is  the  time  required  to  decrease  the  surface  temperature  from  720®F  to  80°F.  For 
APC-2  the  temperatures  were  also  recorded  at  the  midpoint  inside  the  composite. 

After  having  been  cooled  to  room  temperature  the  crystallinities  of  the  specimens 
were  determined  by  differential  scanning  calorimetry  (DSC).  A  small  sample  taken  from  the 
specimen  (~40  mg)  was  enclosed  in  an  aluminum  saunple  pan  and  placed  into  the  DSC  cell. 
The  heat  flow  in  the  cell  as  a  function  of  temperature  was  measured  at  a  constant  heating 
rate  of  36®F/min.  The  crystallinity  was  determined  from  the  expression  (Figure  26) 

^  (52) 

tlu 

where  Hf  vs  the  net  amount  of  heat  absorbed  on  heating  from  the  glass  transition  temper¬ 
ature  (293°F)  to  720° F  as  measured  by  the  DSC.  is  the  ultimate  heat  of  fusion  and  is 
130  J/gram  [17].  The  crystallinities  of  some  of  the  PEEK  150P  specimens  were  also  deter¬ 
mined  by  wide  angle  X-ray  diffraction  [15].  The  crystallinities  measured  by  DSC  and  X-ray 
agreed  closely.  In  evaluating  the  crystallinities  of  the  APC-2  composites  adjustments  were 
made  for  the  fiber  content  using  the  expression 

c  =  I(c)„,„„,rfl/(1VP)  (53) 

where  WP  is  the  weight  percent  of  the  polymer  in  the  composite.  Since  the  polymer  density 
depends  on  the  crystallinity,  an  iterative  procedure  was  used  to  calculate  WP  [15].  This 
procedure  was  based  on  a  nominal  fiber  volume  fraction  of  0.59. 

The  crystallinities  of  1/8  inch  thick  PEEK  150P  plates  specimens  manufactured  at 
three  different  cooling  rates  are  presented  in  Figure  27.  The  temperatures  recorded  at  the 
middle  of  3/16  inch  thick  APC-2  composites  are  shown  in  Figure  28. 

The  crystallinities  and  the  center  temperatures  were  also  computed  by  the  crystallinity 
submodel.  The  material  properties  used  in  the  calculations  are  listed  in  the  Appendices. 
The  calculated  values  of  the  temperatures  and  crystallinities  and  temperatures  are  included 
in  Figures  27  and  28.  There  are  good  agreements  between  the  data  and  the  results  of  the 
model. 
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TIME  ( TEMPERATURE) 


Figure  26. 


Illustration  of  the  results  obtained  on  heating  a  PEEK  150P  polymer  or  APC-2 
composite  sample  in  the  DSC. 
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CRYSTALLINITY,  wt.  %  TEMPERATURE  (°F ) 


TIME  (min) 


COOLING  RATE 


Figure  27.  Average  crystallinities  of  PEEK  150P  plates  cooled  at  different  rates.  Com¬ 
parisons  between  the  measured  and  calculated  crystallinities.  Top:  Cooling 
cycle  used  in  the  manufacture  of  the  plates. 
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CENTER  TEMPERATURE  (°F)  CENTER  TEMPERATURE  (°F  ) 


TIME,  t  (sec) 


Figure  28.  Temperature  in  the  middle  of  APC-2  plates  cooled  (heated)  at  different  rates. 
Comparison  between  the  data  and  the  results  of  the  model. 
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Section  V 


NUMERICAL  PROCEDURE 

Solution  to  the  impregnation,  intimate  contact,  autohesion,  and  crystallinity  submodels 
may  be  obtained  individually  by  the  methods  described  in  conjunction  with  each  submodel. 
Some  of  the  solutions  are  simple,  while  other  require  extensive  numerical  computations.  For 
convenience,  all  of  the  four  submodels  were  incorporated  in  a  user  friendly  computer  code 
designated  as  PLASTIC.  The  input  parameters  required  by  this  code  are  summarized  in 
Tables  1  through  3.  The  output  provided  by  the  code  is  given  in  Table  4  and  in  Figure  29. 
With  this  code  each  processing  step  can  be  studied  individually.  The  code  can  also  be  used 
to  analyze  processes  where  two  or  more  phenomena  occur  simultaneously.  The  model,  and 
the  corresponding  code,  should  thus  be  useful  in  the  study  of  the  processing  of  thermoplastic 
matrix  composites,  and  in  establishing  the  processing  parameters  (temperature,  pressure, 
time)  most  suitable  in  a  given  application. 
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Table  4 


OUTPUT  PROVIDED  BY  THE  THERMOPLASTIC  PROCESSING  MODEL 

AND  THE  PLASTIC  CODE 

Impregnation  Submodel 

1)  Degree  of  impregnation  as  a  function  of  time. 

2)  Time  required  to  complete  the  impregnation. 

Consolidation  Submodel 

3)  Degree  of  intimate  contact  as  a  function  of  time. 

4)  Time  required  for  complete  intimate  contact. 

5)  Degree  of  autohesion  as  a  function  of  time. 

6)  Time  required  for  complete  autohesion. 

7)  Degree  of  bonding  as  a  function  of  time. 

8)  Time  required  for  complete  bonding. 

Crystallinity  Submodel 

9)  Temperature  inside  the  material  as  a  function  of  position  and  time  during  cooling. 

10)  Crystallinity  inside  the  material  as  a  function  of  position  and  time  during  cooling. 

11)  Crystalhnity  as  a  function  of  position  inside  the  material  after  completion  of  the  pro¬ 
cessing. 

12)  Average  crystallinity  of  the  material  after  processing. 

13)  Mechanical  properties  as  a  function  of  position  inside  the  material. 
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Figure  29.  Illustration  of  the  output  of  the  thermoplastic  processing  model. 


Section  VI 


MECHANICAL  PROPERTIES 

The  mechanical  properties  of  thermoplastic  matrix  composites  can  be  related  directly  to 
the  crystalhnity  of  the  polymer.  Thus,  once  the  variations  of  the  mechanical  properties  with 
the  crystalhnity  are  estabhshed  by  tests,  the  mechanical  properties  resulting  from  a  given 
cooling  rate  can  simply  be  determined  from  the  model,  as  follows.  First  the  crystalhnity 
is  calculated  by  the  model.  Second,  the  mechanical  property  is  calculated  from  the  known 
crystalhnity — property  relationship.  Mechanical  properties  of  PEEK  150P  polymer  and 
APC-2  composite  as  functions  of  crystallinity  were  measured  during  this  investigations. 
These  tests  and  the  results  are  described  in  the  foUowing  sections  (Sections  VII  and  VIII). 
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Section  VII 


MECHANICAL  PROPERTIES— TEST  METHODS 

The  test  specimens  were  cut  out  of  plates  manufactured  by  the  procedures  described 
below. 

Manufacturing  PEEK  150P 

The  PEEK  150P  plates  were  made  by  placing  the  polymer,  in  powder  form,  in  an 
aluminum  mold.  Two  molds  were  used;  one  was  10  x  7  x  1/2  inch,  and  the  other  was 
6  X  l|  X  0.2  inch  (Figures  30  and  31).  The  first  mold  was  used  to  achieve  crystallinities 
above  25  percent.  The  second  mold  was  used  to  make  plates  with  lower  crystallinity. 

The  10-inch  mold  consisted  of  a  base  and  a  cover,  as  shown  in  Figure  30.  The  PEEK 
150P  powder  was  poured  into  the  mold.  Then,  with  the  cover  removed,  the  mold  was  placed 
into  an  oven  and  heated  to  750®F.  This  temperature  was  chosen  because  local  order  has  been 
reported  to  remain  in  the  PEEK  melt  unless  the  polymer  is  heated  above  its  equilibrium 
melting  temperature  [17].  Blundell  and  Osborn  estimated  this  temperature  to  be  743®F 
[18].  It  is  very  difficult  to  achieve  low  crystallinity  if  the  polymer  is  melted  at  too  low  a 
temperature. 

The  mold  was  kept  at  this  temperature  for  two  hours  to  allow  air  and  vapor  bubbles 
to  escape  from  the  polymer.  The  cover  was  placed  on  top  of  the  polymer  melt  and  the  mold 
was  kept  at  750‘’F  for  an  additional  15  minutes.  The  polymer  was  then  cooled  to  room 
temperature  by  blowing  air  past  the  mold.  The  cooling  rate  was  controlled  by  regulating  the 
air  flow  rate.  During  cooling,  the  temperature  was  monitored  by  a  Type  J  thermocouple 
located  in  the  center  of  the  polymer  plate.  Because  of  the  long  time  at  high  temperature,  the 
top  surface  of  the  PEEK  plate  was  covered  with  a  black  film,  presumably  degraded  PEEK. 
This  film  was  removed  by  sanding  before  further  processing.  The  bottom  surface  was  also 
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HEATING  PLATE 


HEATING  PLATE 
8  in  X  8  In  X 1 1n 


Figure  31.  Mold  used  for  preparing  specimens  of  less  than  20  percent  crystallinity. 
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sanded  away  because  it  contained  many  voids.  The  interiors  of  the  plates  were  virtually  free 
of  voids. 

The  6-inch  mold  was  provided  with  heaters  as  well  as  with  built-in  cooling  channels 
(Figure  31).  Details  of  this  apparatus  were  described  elsewhere  [19]  and  are  not  repeated 
here.  Only  the  essential  features  of  the  mold  are  given  below.  Six  7/16  inch  diameter  holes 
were  provided  in  the  top  and  the  base  of  the  mold.  Cold  water  could  be  passed  through 
these  holes  to  provide  rapid  cooling.  The  cooling  rates  were  adjusted  by  regulating  the 
water  flow  up  to  30  gallons  per  minute.  Heating  plates,  containing  1500-watt  cartridge 
heaters,  were  placed  above  and  below  the  mold.  The  heating  plate-mold  assembly  could  be 
placed  into  a  press  which  could  exert  a  pressure  on  the  material  inside  the  mold — polymer 
or  composite — during  the  entire  manufacturing  process,  i.e.,  during  heating  and  cooling. 

Prior  to  processing  in  the  6-inch  mold,  the  PEEK  150P  powder  was  dried  in  an  oven 
at  300° F  for  three  hours.  The  mold  was  then  filled  with  the  powder  to  a  depth  of  about 
0.15  inch.  The  cover  was  placed  on  the  mold  and  the  mold  was  inserted  into  the  press. 
The  temperature  of  the  mold  was  raised  to  390°F  under  a  pressure  of  1500  psig.  After  5 
minutes,  PEEK  150P  powder  was  added  to  increase  the  thickness  of  the  plate.  With  the 
cover  in  place  the  mold  was  remounted  in  the  press.  A  pressure  of  1500  psig  at  390°F  was 
again  applied  for  5  minutes.  Then  the  pressure  was  removed  and  the  mold  temperature  was 
raised  to  750°F  and  held  for  20  minutes  without  pressure  being  applied.  The  data  of  Lee 
and  Porter  [20]  and  Ma  et  al.  [21]  suggest  that  the  time  at  high  temperature  was  not  long 
enough  for  significant  degradation  to  occur  in  the  polymer.  At  the  end  of  20  minutes  the 
heaters  were  turned  off,  and  insulating  plates  were  inserted  between  the  heating  plates  and 
the  mold.  The  pressure  was  raised  to  300  psig  and  the  cooling  water  flow  was  started.  The 
mold  temperature  dropped  to  room  temperature  in  5  to  10  seconds. 

During  the  entire  manufacturing  process  the  temperature  was  monitored  by  an  Omega 
Type  J  subminiature  thermocouple  mounted  on  the  inside  center  of  the  base  of  the  mold. 

Manufacturing  APC-2 

The  APC-2  composite  plates  used  in  this  study  were  provided  by  ICI  and  the  Swedish 
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National  Defense  Institute.  These  plates  were  0.05  and  0.2  inch  thick  (8  and  34  plies)  and 
contained  unidirectional  AS4  graphite  fibers.  Some  of  the  test  specimens  were  cut  directly 
from  these  plates.  In  addition,  sections  of  the  plates  were  reprocessed  in  the  6  inch  mold 
described  above  to  obtain  test  specimens  with  lower  levels  of  crystallinity.  The  procedure 
was  essentailly  the  same  as  that  used  for  the  PEEK  150P  polymer.  The  mold  was  placed 
into  the  press,  a  pressure  of  200  psig  was  applied,  and  the  temperature  was  raised  to  750°F 
in  20  to  25  minutes.  The  temperature  was  kept  at  750°F  and  the  pressure  at  200  psig  for 
10  minutes.  Again,  this  short  time  at  high  temperature  is  not  expected  to  allow  the  matrix 
to  degrade  significantly.  The  heaters  and  the  pressure  were  then  turned  off  and  insulating 
plates  were  inserted  between  the  heating  plates  and  the  mold;  the  pressure  (200  psig)  was 
reapplied,  the  water  was  turned  on,  and  the  plates  were  cooled  to  room  temperature  within 
3  to  10  seconds. 

Annealing  PEEK  150P  and  APC-2 

Some  of  the  PEEK  150P  polymer  and  APC-2  composite  plates  were  annealed  as  follows. 
Plates  were  processed  in  the  6-inch  mold  using  a  high  cooling  rate,  resulting  in  plates  with 
low  (below  20  percent)  crystallinities.  After  the  plate  was  cooled  to  room  temperature  the 
mold  was  reheated  in  20  to  25  minutes  to  480®F  [18].  During  this  reheat  no  pressure  was 
applied  for  the  PEEK  150P  polymer  and  200  psig  wais  applied  for  the  composite.  Once 
480°F  was  reached,  the  temperature  was  kept  at  this  level  for  60  minutes.  The  mold  was 
then  cooled  to  room  temperature. 

Differential  Scanning  Calorimetry  and  X-ray  Inspection 

The  crystallinities  of  polymer  and  composite  samples  were  determined  using  differen¬ 
tial  scanning  calorimetry  (DSC)  and  the  peak  area  calculation  method  [22].  Although  for 
PEEK  of  low  crystallinity  the  density  gradient  technique  is  more  accurate  than  DSC  [23], 
the  presence  of  voids  in  composites  makes  density  measurements  impractical  for  APC-2. 
Therefore,  DSC  was  used  in  this  study. 

A  small  sample  taken  from  the  plate  (~10  mg)  was  enclosed  in  an  aluminum  sample 
pan  and  placed  into  the  DSC  cell.  The  heat  flow  in  the  cell  as  a  function  of  temperature 


59 


was  measured  at  a  constant  heating  rate  of  36°F/minute.  The  crystallinity  was  determined 


from  the  expression 


c  = 


Ht 

H\]lt 


(54) 


where  Ht  is  the  net  amount  of  heat  absorbed  on  heating  from  the  glass  transition  temper¬ 
ature  (293°F)  to  720°F  as  measured  by  the  DSC.  i^ULT,  the  ultimate  heat  of  crystallization 
of  100  percent  crystalline  polymer,  was  taken  to  be  130  J/g  [18]. 

In  evaluating  the  crystallinities  of  the  APC-2  composite,  adjustments  were  made  for 
the  fiber  content  using  the  expression 

,dl  (55) 


=  [(c) 


measurer 


where  WP  is  the  weight  percent  of  the  polymer  in  the  composite.  Since  the  density  of 
the  polymer  depends  on  the  crystalUnity,  an  iterative  procedure  was  used  to  calculate  the 
weight  percent  polymer  in  the  composite.  This  procedure  was  based  on  a  nominal  fiber 
volume  fraction  of  0.59. 

The  measured  crystallinities  versus  cooling  rates  are  given  in  the  next  section  (Results). 
It  is  noted  here  that  the  crystallinities  as  a  function  of  cooling  rate  of  PEEK  150P  measured 
in  this  study  by  DSC  agree  reasonably  well  with  those  obtained  by  Velisaris  and  Seferies 
[24]  for  PEEK  450P  using  the  density  gradient  technique  (Figure  32).  The  differences  in 
crystallinities  at  high  cooling  rates  may  be  due  to  differences  in  the  material  tested  or 
differences  in  the  measuring  techniques. 

In  addition  to  differential  scanning  calorimetry,  the  crystallinities  of  some  PEEK  150P 
plates  were  determined  by  wide-angle  X-ray  diffraction  using  a  Picker  diffractometer  with 
copper  Ka  radiation.  The  scans  were  analyzed  using  a  method  discussed  by  Young  [25]  and 
applied  to  PEEK  by  Blundell  and  Osborn  [18].  A  straight  line  was  drawn  on  the  diffraction 
scan  between  26  values  of  10°  and  36°  to  subtract  background  radiation.  {6  is  the  angle 
of  incidence  of  the  incoming  X-ray  beam.)  The  area  under  the  diffraction  curve  and  above 
this  line  was  then  divided  into  the  crystalline  and  amorphous  peaks,  as  illustrated  in  Figure 
35.  The  ratio  of  the  area  of  the  crystalline  peak  to  the  area  of  the  amorphous  peak  equals 
the  weight  ratio  of  crystalline  to  amorphous  material.  The  crystallinities  measured  by  X-ray 
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Figure  32.  Crystallinity  of  PEEK  as  a  function  of  cooling  rate.  Data  of  Velisaris  and  Seferis 
[24]  for  PEEK  450P  and  data  from  the  present  work  for  PEEK  150P. 
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X-Ray  Intensity 


Figure  33. 


10  36 


Diffraction  Angle  20  (degrees) 


A  typical  X-ray  scan  of  semicrystalline  PEEK  150P,  showing  the  background 
radiation  and  the  scattering  from  crystalhne  and  from  amorphous  domains  in 
the  sample. 
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diffraction  and  by  DSC  agree  closely,  both  for  specimens  processed  by  cooling  only  and  for 
annealed  specimens  (Figure  34). 

It  is  recognized  that  the  crystallinity  may  vary  across  the  plate.  In  order  to  estimate  the 
nonuniformities  in  crystallinities,  the  crystallinities  across  a  0.2  inch  thick  APC-2  plate  were 
calculated  for  “low”  cooling  rates,  using  the  method  described  in  Reference  26.  The  results, 
given  in  Figure  35,  show  that  the  crystallinities  are  very  uniform  across  the  plates.  Hence, 
for  low  cooling  rates  nonuniformities  in  crystallinity  are  not  expected  to  affect  significantly 
the  properties  measured  in  this  investigation.  At  “fast”  cooling  rates  there  might  be  a 
crystallinity  gradient  across  the  plate.  For  example,  DSC  measurements  on  0.06  inch  thick 
PEEK  150P  plates,  cooled  at  ll,000®F/min,  showed  that  the  crystallinity  varied  from  16 
percent  at  the  surface  to  23  percent  in  the  center.  Similar  variations  were  observed  by 
Velisaris  and  Seferis  [24]  in  0.15  inch  thick  PEEK  450P  plates  cooled  at  50°F/min. 
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Crystallinity,  c 
by  X-Ray  (percent) 


Figure  34. 


(percent) 


A  comparison  of  the  crystallinities  of  PEEK  150P  as  measured  by  DSC  and  by 
X  ray. 
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Figure  35.  Calculated  distributions  of  crystallinity  in  thick  APC-2  plates  at  different  cooling  rates 


Section  VIII 


MECHANICAL  PROPERTIES— RESULTS 

Tests  were  performed  to  relate  the  mechanical  properties  of  PEEK  150P  polymer  and 
APC-2  composite  to  the  crystallinity  of  the  polymer.  The  results  of  the  tests  are  summarized 
below.  It  is  emphasized  that  these  results  are  meant  primarily  to  aid  in  engineering  analyses. 
No  attempt  has  yet  been  made  to  generate  data  of  sufficient  accuracy  for  a  study  of  the 
detailed  molecular  processes  and  morphologies  involved. 

PEEK  150P  Polymer 

The  PEEK  150P  polymer  test  specimens  were  cut  with  a  band  saw  from  plates  made 
by  the  procedures  described  in  the  previous  section.  Six  types  of  mechanical  properties  were 
measured  at  different  crystallinities:  tensile,  compressive,  and  shear  strengths,  tensile  and 
shear  moduli,  and  mode  I  fracture  toughness. 

The  tensile  properties  were  measured  with  6-in  long  (1/4  inch  by  1/10  inch)  dog-bone 
specimens  (Figure  36).  Data  were  taken  at  a  crosshead  speed  of  0.2  inch/minute.  The 
moduli  were  determined  from  strains  measured  with  a  1  inch  gauge  length  extensometer. 
The  results  are  presented  in  Table  6  and  Figure  37,  where  each  point  is  the  average  of  at 
least  four  tests.  The  tensile  strengths  St  and  tensile  moduli  Et  of  specimens  prepared  by 
coohng  only  can  be  fitted  by  the  expressions 

5*  =  8.2 -k  0.17  c  (56) 

Et  =  404  -1-  4.0  c  +  0.075  (57) 

where  c  is  the  crystallinity  in  percent.  Jones  et  al.  [27]  reported  a  yield  strength  of  15.1  ksi 
for  PEEK  of  unspecified  grade  and  crystallinity.  This  value  is  within  the  range  of  the  tensile 
strengths  found  in  the  current  work  (Table  6). 

The  shear  properties  of  PEEK  150P  were  determined  by  the  losipescu  method  [35]. 
The  geometry  of  the  test  specimens  is  shown  in  Figure  38.  The  tests  were  performed  at  a 
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Table  5 

GEOMETRIES  OF  THE  TEST  SPECIMENS 


Specimen 

Length 

L  (in) 

Width 
w  (in) 

Thickness 
h  (in) 

Notch  Length 
(in) 

Remarks 

Tensile 

6 

1/4 

1/10 

- 

PEEK150P 

Shear 

3 

3/4 

1/10 

0.15 

PEEK150P 

Compressive 

3/8 

3/8 

«1/10 

- 

PEEK150P 

Single  edge 
notched  fracture 

5 

0.06-0.15 

3/4 

PEEK150P 

Edge  notched 
flexure 

4 

3/4 

0.2 

1 

APC-2 

Center  notched 

6 

0.05 

«l/8 

APC-2 

Center  notched 
(off  axis) 

6 

U 

0.05 

a  sin  ^  j 

APC-2 

Rail  shear 

6 

6 

0.05 

1/2 

APC-2 
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MEASURED  PROPERTIES  OF  PEEK  150P 
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Tensile  Modulus 


Crystallinity,  c  (percent) 

Figure  37.  Tensile  properties  of  PEEK  150P  as  functions  of  crystallinity.  The  error  bars 
represent  the  standard  deviation  of  the  data. 
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loading  rate  of  0.24  inch/minute.  The  moduli  were  evaluated  using  strain  gauges  attached 
to  the  specimens.  The  results  are  given  in  Table  6  and  Figure  38,  where  each  point  is  the 
average  of  four  measurements.  The  shear  strength  S',  and  modulus  Eg  data  of  specimens 
prepared  by  coohng  only  were  fitted  by  the  expressions 

S,  =  4.9  + 0.10  c  + 0.0015  (58) 

Es  =  140  +  1.59  c  (59) 

where  again,  c  is  in  percent.  The  foregoing  data  show  that  the  strength  and  stiffness  in  both 
tension  and  shear  increase  with  crystallinity.  The  likeliest  reason  for  the  rise  in  modulus  is 
the  high  stiffness  of  the  crystal  lattice  [29]. 

The  strength  of  PEEK  150P  in  compression  was  measured  using  3/8  inch  square  by 
approximately  1/10  inch  thick  specimens  (Figure  36).  The  specimens  were  compressed  be¬ 
tween  two  plates  at  a  rate  of  0.2  inch/minute.  The  general  form  of  the  load  versus  deflection 
curves  is  shown  in  Figure  39.  As  this  figure  illustrates,  it  is  difficult  to  define  a  unique  com¬ 
pression  strength.  Therefore,  the  compression  strength  was  defined  as  the  point  where  the 
slopes  of  the  first  two  adjacent  segments  intersect.  The  compression  strengths  thus  obtained 
are  presented  in  Table  6  and  Figure  40.  In  this  figure  and  table  each  point  is  the  average  of 
at  least  five  data.  The  compression  strength  Sc  data  in  Figure  39  for  specimens  prepared  by 
coohng  only  can  be  represented  by  the  expression 

=  18  +  0.22  c  (60) 


where  c  is  in  percent. 

It  is  seen  that  the  strength  in  compression  changes  very  little  with  crystallinity  for 
specimens  prepared  by  cooling  only.  However,  due  to  the  highly  nonlinear  behavior  of  the 
material  during  compression  (Figure  39),  these  data  must  be  interpreted  with  caution.  In 
fact,  the  shape  of  the  load- deflection  curves  suggests  that  some  change  in  crystal  structure 
occurs  at  the  point  where  the  curve  bends  over  sharply.  X-ray  scans  of  compressed  specimens 
did  indeed  reveal  a  considerable  reduction  in  crystallinity  after  the  specimen  was  loaded  sig¬ 
nificantly  past  the  deflection  point.  Specimens  of  16  percent  crystallinity  before  compression 


71 


Shear  Strength  Shear  Modulus 


Figure  38. 
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Shear  properties  of  PEEK  150P  as  functions  of  crystallinity.  The  error  bars 
represent  the  standard  deviation  of  the  data. 
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Load 


Deflection 


Figure  39.  A  typical  load-deflection  curve  during  compression  tests  of  PEEK  150P. 
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Compression  Strength 


Crystallinity,  c  (percent) 


Figure  40.  Strength  in  compression  of  PEEK  150P  a  function  of  crystallinity.  The  error 
bars  represent  the  standard  deviation  of  the  data. 
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were  less  than  5  percent  crystalline  thereafter,  while  specimens  at  27  percent  crystalhnity 
before  compression  were  afterwards  only  16  percent  crystalline. 

Finally,  the  mode  I  fracture  toughness  of  PEEK  150P  at  various  degrees  of  crystallinity 
was  measured  using  single  edge  notched  specimens  (Figure  36)  pulled  to  failure  at  0.2  inch 
per  minute.  A  fracture  toughness  K  is  calculated  by  [30] 

PY^ 
hw 

where  P  is  the  failure  load  in  Ibf,  a  is  the  initial  crack  length  in  inches,  h  is  the  specimen 
thickness  in  inches,  and  w  is  the  specimen  width  in  inches.  Y  is  the  dimensionless  finite 
width  correction  factor  computed  by 

Y  =  1.99  -  0.41  4- 18.7  -  38.48  +  53.85  (62) 

If  the  test  specimen  is  in  plane  strain,  the  measured  fracture  toughness  K  is  the  actual 
fracture  toughness  Kjc]  if  plane  strain  conditions  are  not  met  during  the  test  then  K  is  an 
apparent  fracture  toughness  symbolized  by  Kq. 

The  plane  strain  condition  is  met  if  the  ratio  R  is  less  than  unity  [30] 

R  = - 7=  63 

5.^/A 

As  can  be  seen  from  Table  7,  only  specimens  of  33  and  39  percent  crystallinity  were 
in  plane  strain;  therefore,  only  at  these  crystallinities  are  the  measured  toughnesses  for 
PEEK  150P  equal  to  Kjc-  At  lower  crystallinities  the  measured  fracture  toughness  is  an 
apparent  fracture  toughness  Kq.  A  sufficient  increase  in  specimen  thickness  would  produce 
plane  strain  in  these  specimens.  Unfortunately,  it  is  extremely  difficult,  if  not  impossible,  to 
make  thick  specimens  of  low  and  uniform  crystalhnity.  The  measured  values  of  the  fracture 
toughnesses  are  given  in  Table  6  and  Figure  41,  where  each  point  is  the  average  of  at  least 
five  data.  The  decrease  in  toughness  with  increasing  crystallinity  is  quite  pronounced.  A 
similar  trend  has  been  observed  for  isotactic  polystyrene  [31]. 

The  fracture  toughness  data  obtained  in  this  study  may  be  compared  with  the  data  of 
Karger-Kocsis  and  Friedrich  [32],  Klei  [33],  and  Jones  et  al.  [27].  Karger-Kocsis  and  Friedrich 
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Table  7 

THE  PLANE  STRAIN  CRITERION 
FOR  FRACTURE  TOUGHNESS  TESTING  OF  PEEK  150P 


Crystallinity 
c  (percent) 

Measured 

Fracture  Toughness 
K  (ksiym) 

Tensile 

Strength 

Si{ksi) 

Thickness 
h  (inches) 

R 

Remarks 

16 

10.4 

10.5 

0.06 

6.4 

II 

26 

5.4 

13.4 

0.13 

1.8 

K  =  Kq 

27“ 

5.6 

16.5 

0.06 

2.2 

k  =  Kq 

33 

3.3 

14.1 

0.15 

0.96 

K  =  Kic 

39 

2.6 

14.7 

0.15 

0.72 

o 

II 

“Annealed. 
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Mode  I  Fracture  T oughness 
K  (ksi  Vin) 


Figure  41. 


Crystallinity,  c  (percent) 


Measured  Mode  I  fracture  toughness  of  PEEK  150P  as  a  function  of  crystallinity. 
The  error  bars  represent  the  standard  deviation  of  the  data. 


77 


measured  the  fracture  toughness  of  PEEK  450P  of  24  and  30  percent  crystallinities.  At  the 
lower  crystallinity  the  fracture  toughnesses  were  6  to  6.3  ksi’>/m  (at  0.04  inch/min  strain  rate) 
and  0.9  to  2.1  ksiv^  (at  0.4  inch/min).  At  the  higher  crystallinity  the  fracture  toughnesses 
were  5.2  and  0.8  ksiVm  at  0.04  and  0.4  inch/min  strain  rates,  respectively.  These  fracture 
toughness  values  are  similar  to  those  obtained  in  the  present  investigation  (see  Table  6  and 
Figure  41). 

Klei  [33]  reported  the  fracture  toughnesses  of  PEEK  at  21  to  25  percent  crystallinity. 
The  grade  of  PEEK  used  in  the  tests  was  unspecified.  The  toughnesses  measured  by  Klei 
(11  to  15  ksiv^)  are  somewhat  higher  than  those  presented  here  (2.5  to  11  ksi-^/in).  The 
differences  in  toughness  obtained  by  Klei  and  by  us  may  be  due  to  differences  in  material, 
in  testing  speeds  (0.02  inch/min  versus  0.2  inch/min),  or  sample  thickness  (0.1  inch  versus 
0.15  inch). 

Jones  et  al.  [27]  did  not  specify  either  the  grade  of  PEEK  or  the  crystallinity  used 
in  their  tests.  They  measured  the  toughness  of  PEEK  as  4.1  ksi-y/m  (at  a  testing  speed  of 
3.9  inch/min),  a  value  comparable  to  those  presented  here. 

We  measured  the  tensile,  compressive,  and  shear  properties  and  the  fracture  toughness 
for  annealed  as  well  as  for  cooled-only  specimens.  The  results  are  included  in  Table  6  and  in 
Figures  37,  38,  40,  and  41.  From  these  data,  it  appears  that  annealed  specimens  are  stronger 
in  both  tension  and  compression  than  cooled-only  specimens  of  comparable  crystallinity. 
However,  the  shear  strengths,  the  tensile  and  shear  moduli,  and  the  fracture  toughnesses  of 
annealed  specimens  are  similar  to  those  of  cooled-only  specimens  of  the  same  crystallinity. 

APC-2 

For  APC-2  composite,  the  Mode  I  and  Mode  II  fracture  toughnesses  {Kjc  and  Kuc) 
and  fracture  energies  {Gic  and  Guc)  as  well  as  the  mixed  mode  stress  intensity  factors  {Kj 
and  Kji)  were  determined  for  crystalhnities  ranging  from  0  to  33  percent. 

The  Mode  I  fracture  toughnesses  were  measured  using  center  notched  specimens,  with 
the  notch  being  parallel  to  the  fibers  and  normal  to  the  load  (Figure  42).  The  notching 
procedure  is  given  in  [34].  Tensile  load  was  applied  at  a  crosshead  speed  of  0.07  inch/minute. 
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Figure  42.  Geometries  of  APC-2  test  specimens.  The  specimen  dimensions  are  specified  in 
Table  5.  Top  left:  center  notched  (CN),  top  right:  rail  shear  (RS),  bottom: 
edge  notched  fracture. 


(64) 


The  fracture  toughness  was  calculated  by  the  expression  [35] 

PY-s/^ 


KlC  = 


hw 


where  P  is  the  load  at  which  failure  occurs  and,  as  before,  h  and  w  are  the  thickness  and 
the  width  of  the  specimen.  2a  is  the  initial  notch  length.  Y  is  the  finite  width  correction 
factor  defined  as 

2  3 

r  =  1.77-1-0.227  -0.510  -1-2.7^^^  (65) 

From  the  known  values  of  Kjc,  the  Mode  I  critical  strain  energy  release  rates  Gjc  (re¬ 
ferred  to  here  in  short  as  the  Mode  I  “fracture  energies”)  were  calculated  by  the  formula  [36] 


Gic  =  Kic 


{2ElEtYI'^ 


.Et) 


1/2 


l^LT  + 


El 

2Glt 


1/2 


(66) 


where  El  and  E'j'  are  the  longitudinal  and  transverse  elastic  ply  moduli,  Glt  i®  the  shear 
ply  modulus  and  i/lt  i®  the  ply  Poisson  ratio.  In  the  calculations  El  and  i^lt  were  taken  to 
be  constant  [37]  and  having  the  values  given  in  Table  8.  The  transverse  and  shear  moduli 
and  the  Poisson  ratio  were  estimated  from  the  measured  polymer  properties  given  in  Table  8 
and  the  micromechanics  formulae  given  in  Table  9.  The  properties  thus  estimated  are  given 
in  Table  10. 


The  Mode  II  fracture  toughnesses  were  measured  using  either  rail  shear  (RS)  or  edge 
notched  fracture  (ENF)  specimens  (Figure  42).  The  rail  shear  test  was  performed  on  spec¬ 
imens  of  high  crystallinity.  Edge  notched  fracture  tests  were  used  with  specimens  at  low 
crystallinities  because  these  specimens  (processed  in  the  6-inch  mold)  were  too  small  for  the 
rail  shear  test. 


For  the  rail  shear  specimens,  load  was  applied  at  a  rate  of  0.03  inch/minute.  The 
fracture  toughness  was  calculated  from  the  measured  failure  load  P  by  [39] 

KlIC  =  -^y/^  (67) 

where  h  and  L  are  the  specimen  thickness  and  length,  and  2a  is  the  initial  notch  length. 
Correction  for  finite  width  was  not  applied  because  such  a  correction  was  shown  to  be 
unnecessary  for  the  type  of  specimens  used  in  this  study  [39]. 
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Table  8 

PROPERTIES  USED  IN  CALCULATING 
THE  ELASTIC  PROPERTIES  OF  APC-2  GIVEN  IN  TABLE  10 


Fiber 

Volume  fraction,  Vf  =  0.58 
Poisson’s  ratio,  i/f  =  0.27 
Transverse  modulus.  Erf  =  2.32  Msi 
Longitudinal  shear  modulus,  Gf  =  1.2  Msi 

Matrix 

Volume  fraction,  Vm  =  0.42 
Matrix  modulus, 

Stress  Partitioning  Parameters:  tjlt  —  0-2i  VT  =  0.5. 
Composite 

Poisson’s  ratio,  i/it  —  0-28  ' 

Longitudinal  modulus,  Ei  =  19.4  Msi 


measured 
from  [37] 
from  [37] 
from  [37] 


measured 
from  Table  6 
from  procedure  in  [37] 


from  [37] 
from  [37] 


Table  9 

MICROMECHANIC  EQUATIONS  USED 
IN  CALCULATING  ELASTIC  PROPERTIES  OF  APC-2 

(Tsai-Hahn  [38]) 

Longitudinal  Shear 
Modulus,  Git 

Transverse  Tensile 
Modulus,  Et 

Poisson’s  Ratio  of 
the  Matrix,  i/m 

Shear  Modulus  of 
the  Matrix,  Gm 

E  Tensile  modulus 

G  Shear  modulus 

V  Volume  fraction 

T]  Stress  partitioning  parameter  in  the  Tsai-Hahn  equations 

1/  Poisson’s  ratio 

Subscripts 

/  Fiber 

m  Matrix 

L  Longitudinal  direction 

T  Transverse  direction 


1 _ 


1  _  1  _ 

Et  Vf-^-TJTVm  \ETf 


(JL. 

\ETf 


"*  Vm 


G - Em. 
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Table  10 

ELASTIC  PROPERTIES  USED  IN  CALCULATING 
THE  FRACTURE  ENERGY  G/c 
AND  FRACTURE  TOUGHNESS  Knc  OF  APC-2 


Crystallinity 
c  (percent) 

Transverse  modulus 
Ej’(Afsi) 

Shear  modulus 
Gix{Msi) 

0-20 

1.19 

0.71 

25 

1.21 

0.72 

28 

1.29 

0.74 

33 

1.37 

0.78 

Ex  and  Gix  were  calculated  using  the  equations  in  Table  9 
and  the  data  in  Table  8. 
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Load  was  applied  to  the  HNF  specimens  at  the  rate  of  0.1  inch/miimte.  The  fracture 
energy  Gjic  was  computed  from  the  measured  failure  load  P  by  the  expression  [40] 


9a'^P^C 

~  2w{2L^+5a^) 


(68) 


Again,  L,  w,  and  a  represent  the  length  and  width  of  the  specimen  and  the  initial  notch 
length.  C  is  the  compliance  due  to  bending  and  is  the  ratio  of  the  measured  deflection  at 
the  loading  position  divided  by  the  load.  The  mode  II  fracture  toughness  Kuc  is  related  to 
the  Mode  II  fracture  energy  Gjic  by  [36] 

K„c  =  [Gnc^E,)  ‘''1(1;)  -  "it  +  |  (69) 

The  meanings  of  the  symbols  are  the  same  as  for  Equation  66.  Knc  values  were  calculated 
from  this  equation,  with  the  moduli  Ei,  Et,  Git  and  the  Poisson  ratio  vlt  determined  in 
the  same  manner  as  discussed  above  (see  Equation  66). 


The  Mode  I  and  Mode  II  fracture  toughnesses  and  fracture  energies  as  functions  of 
crystallinity  are  presented  in  Table  11  and  Figures  43  through  44.  Each  data  in  this  table 
and  figures  is  the  average  of  at  least  three  tests. 

The  following  observations  can  be  made  from  the  data  in  Figures  43  and  44.  First,  the 
mode  II  fracture  toughnesses  obtained  by  rail  shear  and  edge  notched  fracture  specimen  tests 
agree  very  closely.  Second,  at  corresponding  crystallinities,  the  fracture  toughnesses  obtained 
from  specimens  which  were  cooled  only  were  nearly  the  same  as  the  fracture  toughnesses 
of  specimens  which  were  first  cooled  rapidly  then  heated  (“annealed  ).  This  observation 
suggests  that  neither  the  fracture  toughness  nor  the  fracture  energy  is  sensitive  to  processing 
history,  but  depends  mainly  on  the  value  of  the  crystallinity.  Third,  the  fracture  toughnesses 
and  fracture  energies  decrease  with  increasing  crystallinity. 

The  data  in  Figures  43  and  44  were  fitted  with  curves,  which  can  be  described  by  the 
equations 


=  4.3  -  0.007  c  -  0.001  ksi\/in 

Knc  =  9-2  -  0.005  c  -  0.002  c^  ksiv^ 
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(70) 

(71) 


Table  11 

FRACTURE  PROPERTIES  OF  APC-2  AS  A  FUNCTION  OF  CRYSTALLINITY 


Fracture  Toughness  Fracture  Energy 

Cooling  rate  Crystallinity  Mode  I  Mode  II  Mode  I  Mode  II 

(°F/niin)  c  (percent)  Kic{ksi\fm)  Knc{ksi\/m)  G/c  (Ibf/in)  G/yc  (Ibf/in) 


18000 

0 

4.3 

9.2 

11.5 

12.1 

6500 

0 

4.4 

8.9 

11.8 

12.7 

4300 

14 

4.3 

9.0 

11.3 

12.3 

1700  to  430°F 

90  from  430*’F 

20 

4.1 

8.4 

10.2 

10.7 

18000+  annealing 
for  60  min. 
at  480°F 

25 

3.3 

8.3 

6.7 

10.3 

72 

28 

3.4 

8.1 

6.8 

9.7 

27 

33 

3.2 

6.8 

5.6 

6.7 
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Crystallinity,  c  (percent) 

Figure  43.  Mode  I  and  Mode  II  fracture  toughness  of  APC-2  as  a  function  of  crystallinity. 

The  error  bars  represent  the  standard  deviation  of  the  data.  CN  =  center 
notched,  ENF  =  edge  notched  fracture,  RS  =  rail  shear  specimen  (see  Figure 
421. 


Fracture  Energy  Fracture  Energy 

G||c  (Ibf/in)  ®lc 


0  10  20  30  40 

Crystallinity,  c  (percent) 


Figure  44.  Mode  I  and  Mode  II  fracture  energy  of  APC-2  cis  a  function  of  crystallinity.  The 
error  bars  represent  the  standard  deviation  of  the  data.  CN  =  center  notched, 
ENF  =  edge  notched  fracture,  RS  =  rail  shear  specimen  (see  Figure  42). 
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Gic  =  12.01  -  0.022  c  -  0.006 

Ibf/in 

(72) 

Giic  =  12.78  -  0.024  c  -  0.004  c^ 

Ibf/in 

(73) 

where  c  is  the  degree  of  crystallinity  in  percent. 


The  Mode  I  fracture  toughness  obtained  in  this  investigation  can  be  compared  to  the 
value  given  by  Donaldson  [41],  who  measured  the  fracture  toughness  of  APC-1.  Donaldson’s 
result  of  3.4  ksiv^  is  the  same  as  our  valve  obtained  at  28  percent  crystallinity  (see  Table 
11). 


Mode  I  and  Mode  II  fracture  energies  of  APC-2  have  previously  been  reported  in 
References  42  through  45.  Russel  and  Street  [42],  Smiley  and  Pipes  [43]  and  Leach  et  al.  [44] 
obtained  Gjc  values  of  8.8,  9.1,  and  13.8  to  16.5  Ibf/in,  respectively.  These  values  are  similar 
to  those  obtained  during  the  course  of  this  study  (5.6  to  11.8  Ibf/in,  Table  11).  Russell  and 
Street  [42]  and  Carlsson  et  al.  [45]  reported  Gjjc  values  of  10.1  and  10.6  Ibf/in,  respectively. 
These  values  are  also  in  the  range  of  those  values  measured  here  (6.7  to  12.7  Ibf/in,  Table  11). 


Finally,  the  stress  intensity  factors  Kj  and  Ku  were  measured  under  mixed  mode 
loading  conditions,  using  specimens  with  crystallinities  ranging  from  0  to  33  percent.  The 
test  specimen  is  shown  in  Figure  45.  Notches  of  length  2a  were  inserted  parallel  to  the 
fibers.  Specimens  with  fiber  and  notch  angles  of  0  =  45,  30,  and  20  degrees  were  tested  at  a 
crosshead  speed  of  0.07  inch/minute.  The  stress  intensity  factors  at  the  point  of  failure  were 
calculated  from  the  measured  failure  loads  using  the  expressions  [36] 


Ki  = 


y/ifa  sin^  0 


(w)(h) 
p 

—  7 — T77T  sin  6  cos  6 


(74) 

(75) 


{w){h) 

No  finite  width  correction  factors  were  used,  since  none  are  known  to  us  for  off-axis  center- 
notched  specimens.  The  data  are  presented  in  Figure  46.  Each  point  in  this  figure  is  the 
average  of  at  least  three  measurements. 


It  has  been  suggested  in  the  past  that  the  mixed  mode  stress  intensity  factors  for 


composites  can  be  related  by  the  equation  [46] 


\KicJ  \KiicJ 


(76) 
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In  Figure  46  the  above  equation  was  also  included.  As  can  be  seen,  the  quadratic  failure 
criterion  represented  by  Equation  76  describes  the  data  well  over  the  entire  range  of  crys¬ 
tallinities  tested.  This  supports  the  observations  made  by  previous  investigators  [34,  47] 
that  the  mixed  mode  fracture  of  thermoplastic  composites  is  controlled  by  the  total  fracture 
energy. 
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p 


p 


Figure  45. 


Geometry  of  mixed-mode  fracture  specimens  of  APC-2.  The  specimen  dimen¬ 
sions  are  specified  in  Table  5. 
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Stress  Intensity  Factor,  K|  (ksi  Vin) 


Figure  46. 


Stress  Intensity  Factor,  K||  (ksi  Vin) 

Mixed-mode  fracture  toughness  of  APC-2  as  a  function  of  crystallinity.  The 
error  bars  represent  the  standard  deviation  of  the  data. 
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APPENDIX  A 


Material  Properties  Used  in  the  Thermoplastic  Processing  Model 

In  this  appendix  the  material  properties  pertaining  to  the  processing  of  PEEK  150P 
polymer,  T300/PEEK  150P  and  APC-2  composites  are  listed.  Additional  information  on 
the  crystallinity  of  PEEK  150P  polymer  is  provided  in  Appendix  B. 


A-1 


Table  A-1 


VALUES  OF  THE  PARAMETERS  IN  THE  IMPREGNATION 
CONSOLIDATION  AND  CRYSTALLINITY  SUBMODELS 


Impregnation — T300  Fiber 


Initial  tow  radius 
Fiber  diameter 

Average  distance  between  two  fibers 
(channel  height) 

Average  channel  length 

Permeability  of  the  tow 

Surface  tension  (T300/PEEK  150P  at  698°F) 

Viscosity  (PEEK  150P) 


Vo  =  380  fim 
d  =  6  iim 
h  —  2  fim 

L  =  150  fim 
Kp  —  0.7 

<j*  =  0.026  iVm-n=  0.0018  Ibf/ft) 


H  =  1.13  X  10“^° 


19100 


Pa  •  s 


Intimate  Contact 

Geometric  ratio  Wojho  =  1.0 

Geometric  ratio  =  0.3 

Matrix-fiber  viscosity 


/x^/  =  1.14x  10-^2 


26300  ■ 


Pa  •  s 


Autohesion 

Proportionality  constant 


K  =  44.1  exp 


3810 

T{°K) 


,-1/4 


Crystallinity 

PEEK  density 

PEEK  specific  heat 

PEEK  thermal  conductivity 

PEEK  ultimate  heat  of  crystallization 

APC-2  matrix  mass  fraction 

APC-2  ply  thickness 


pm  =  1300  kg/ra  (81.1  Ibm/ft^) 

Cm  =  1340  J/kg/°C  (0.32  BTU/lbm/°F) 
km  =  0.25  J/s/m/°C  (0.145  BTU/hr/ft/°F) 
Hu  =  130  J/g  (56  BTU/lbm) 
rrim  =  0.31 

So  =  0.125  mm  (0.005  in) 
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Table  A- 2 


MECHANICAL  PROPERTIES  OF  PEEK  150P  POLYMER  AND 
UNIDIRECTIONAL  APC-2  COMPOSITE  AS  FUNCTIONS  OF  CRYSTALLINITY  [15] 

PEEK  150P 

Tensile  strength:  St  =  8.2  +  0.17  c  (ksi) 

Tensile  modulus:  Et  =  404  +  4.0  c  +  0.075  (ksi) 

Shear  strength:  Ss  =  4.9  +  0.10  c  +  0.0015  (ksi) 

Shear  modulus:  Es  =  140  +  1.59  c  (ksi) 

Compression  strength:  =  18  +  0.22  c  (ksi) 

Fracture  toughness:  K  =  27.6  *  10”®  ®^^  ^  (ksi)  \/m 

APC-2 

Mode  I  and  Mode  11  fracture  toughnesses  (Kjc  and  Knc)  and  frac¬ 
ture  energies  {Gjc  and  Gjjc) 

Kjc  =  4.3  -  0.007  c  -  0.001  (ksi)Vm 
KilC  =  9-2  —  0.005  c  —  0.002  (ksi)\/m 
Gic  =  12.01  -  0.022  c  -  0.006  (Ibf/in) 

Gnc  =  12.78  -  0.024  c  -  0.004  (Ibf/in) 
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APPENDIX  B 


Crystallinity  of  PEEK  150P  Polymer 


Several  empirical  expressions  have  been  proposed  in  the  past  for  correlating  the  mea¬ 
sured  rate  of  crystallization  with  temperature  [48-54].  We  adopted  the  simple  expression 
proposed  by  Ozawa  [49],  which  is  convenient  and  adequate  for  the  purpose  of  engineering 
analysis.  Ozawa’s  expression  can  be  written  as 


log[-  ln(l  -  Cr)]  =  log  (f>  +  n  log 


or  in  differential  form 


dt  ^  ^''\dTldt)^-^ 

Cr  is  the  relative  crystalhnity  and  is  related  to  the  crystallinity  by 

Hu 


(B.l) 


(B.2) 


(B.3) 


Cr  can  be  determined  from  measurements  performed  in  a  differential  scanning  calorimeter. 
Hu  is  the  theoretical  ultimate  heat  of  crystallization  of  the  polymer.  Ht  is  the  total  heat  of 
crystallization  at  the  given  coohng  rate  (Figure  47).  (/>  is  a  parameter  which  depends  on  the 
temperature  only  and  n  is  a  constant.  The  parameters  (f)  and  n  were  obtained  by  fitting  the 
above  expression  to  the  data  obtained  by  DSC  measurements.  This  procedure  yielded  the 
following  <i>  and  n  values 

(f>  =  exp[— 0.037  T  -f  11.3] 

(BA) 

n  =  0.8 

In  addition,  a  best  fit  to  the  data  yielded  the  following  expression  for  HtIHu 


=  -0.03  In 


+  0.42 


where  dTjdt  is  the  cooling  rate  in  °C/min,  and  Hu  is  [48] 

Hu  =  130  J/g 


(B.S) 


(fl.6) 


Crystallinities  as  a  function  of  cooling  rate  measured  by  different  investigators  for  PEEK 
150P  polymer  and  APC-2  composite  are  given  in  Figure  48.  The  volume  percent  crystallinity 
data  of  Velisaris  and  Seferis  [7,  51]  and  Blundell  et  al.  [52]  were  converted  to  weight  percent, 
using  the  densities  given  in  Reference  [51].  In  Figure  48  crystallinities  computed  by  the 
above  expressions  are  also  shown  (solid  line).  As  can  be  seen  these  expressions  approximate 
the  data  reasonably  well. 
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Figure  47. 
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Illustration  of  the  results  obtained  on  cooling  a  PEEK  150P  pol3Tner  sample  in 
the  DSC. 
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Figure  48.  Crystallinity  of  PEEK  150P  polymer  and  APC-2  composite  as  a  function  of 
cooling  rate  comparisons  between  data  of  Velisaris  and  Seferis  [7,  51],  Blundell 
and  Osborne  [53],  Blundell  et  al.  [52],  Talbott  et  al.  [15],  and  the  model  described 
in  Appendix  B. 
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